Geometric criticality between plaquette phases in integer-spin kagome XXZ 

antiferromagnets 



in 
o 
o 

(N 



in 

(N 



i 

C 

O 
o 



> 
o 

oo 
in 

o 
in 
o 



c 

o 
o 



X 
S3 



Cenke Xu 1 and J. E. Moore 1 ' 2 

1 Department of Physics, University of California, Berkeley, CA 94720 
2 Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720 

(Dated: February 2, 2008) 

The phase diagram of the uniaxially anisotropic s = 1 antiferromagnet on the kagome lattice 
includes a critical line exactly described by the classical three-color model. This line is distinct 
from the standard geometric classical criticality that appears in the classical limit (s — » oo) of the 
2D XY model; the s = 1 geometric T = critical line separates two unconventional plaquette- 
ordered phases that survive to nonzero temperature. The experimentally important correlations at 
finite temperature and the nature of the transitions into these ordered phases are obtained using 
the mapping to the three-color model and a combination of perturbation theory and a variational 
ansatz for the ordered phases. The ordered phases show sixfold symmetry breaking and are similar 
to phases proposed for the honeycomb lattice dimer model and s = 1/2 XXZ model. The same 
mapping and phase transition can be realized also for integer spins s > 2 but then require strong 
on-site anisotropy in the Hamiltonian. 

PACS numbers: 75.10.Jm 75.50 Ee 



I. INTRODUCTION 

The geometric constraints that determine the infinite 
set of ground states in classical highly frustrated magnets 
are known in several celebrated instances to give rise to 
power-law or "critical" correlations at vanishing temper- 
ature. A frustrated magnet is one in which not all in- 
teraction energies can be simultaneously optimized; here 
we follow convention in describing by "highly frustrated" 
the class of frustrated magnets that have an infinite num- 
ber of classical ground states, even after removing global 
symmetries. An example of a frustrated magnet that is 
not highly frustrated is the triangular lattice XY anti- 
ferromagnet, which has a unique ground state once the 
three spins on a single triangle are fixed. 

The classical XY antiferromagnet on the (2D) kagome 
lattice and the classical Heisenberg antiferromagnet on 
the (3D) pyrochlore lattice are special members of the 
highly frustrated class: they show power-law correlations 
determined by a critical two-component height model in 
the first cas o 1 ! 2 ! 3 and by 3D dipolar correlations in the 
second^. A general highly frustrated magnet, such as 
the Ising antiferromagnet on the kagome or triangular 
lattice, has only short-ranged correlations, i.e., an expo- 
nential decay of correlations beyond a finite "correlation 
length." Geometrically critical frustrated magnets can be 
thought of as on the boundary between disordered phases 
with short-ranged correlations and long-ranged ordered 
phases, but rather than being tuned by an external pa- 
rameter, the criticality is a consequence of the fixed ge- 
ometry of the lattice. A major focus of current research 
in frustrated magnetism is how quantum effects in real 
magnets with finite spin S modify these classical critical 



points, leading possibly to spin liquid: 
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emergent 



gauge symmetries*, and novel ordered phases. 

The first goal of this paper is to show that the 
uniaxially anisotropic spin-1 antiferromagnet on the 



kagome lattice has a range of parameters over which 
the ground states are exactly determined by the "three- 
color model"—, which is the same geometrically criti- 
cal model that appears in the classical XY model on 
the same lattic e) 1 ^ 3 . The s = 1 kagome antiferromag- 
net is thus a rare example where exact information can 
be obtained on a frustrated quantum magnet. Uni- 
axial (XXZ) anisotropy is generic in kagome materials 
because their layered quasi-2D structure implies XXZ 
anisotropy. A typical kagome compound is the s = 3/2 
material SrCrg_Ga4+ x Oi^i; nickel kagome compounds 
have s = 1. A recent neutron scattering studj^ of 
the iron jarosite KFe3(OH) 6 (S04)2, which has s = 5/2, 
observed, in addition to expected classical Heisenberg 
kagome physics, several effects induced by Dzyaloshinkii- 
Moriya terms resulting from additional ions that break 
the inversion symmetry of the kagome lattice. 



The second part looks at the plaquette-ordered phases 
separated by this critical line, which serves as a solvable 
starting point for perturbation theory. The critical line 
occurs for reasons quite different from the critical point 
of the classical XY model on this lattice; in particular, 
the s = 1 problem has a nonvanishing gap above the 
ground-state manifold in the thermodynamic limit, un- 
like the large-5 1 case. Section II reviews the connection 
between the classical XY antiferromagnet and the three- 
color model and introduces the XXZ model for spin-1 
antiferromagnets on the kagome lattice. From this we 
derive the existence of a critical line related to the three- 
color model. Section III discusses the ordered phases sep- 
arated by this critical line and the nature of the phase 
transition and concludes by developing a gauge theory 
description in order to compare this system to other re- 
cently studied models of frustrated magnetism. 
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II. THE KAGOME ANTIFERROMAGNET AND 
THE THREE-COLOR MODEL 

The ground states of the classical XY model on the 
kagome lattice (Fig. 1) are equivalen1*ii2i£, up to an over- 
all 0(2) rotation, to the states of the honeycomb lattice 
"three-color model"—. Once one spin is fixed, all other 
spins form angles of or ±2-7r/3 with this spin: labeling 
these three directions with colors, the ground states have 
one spin of each color around each small triangle (Fig. 1). 
This is equivalent to coloring the bonds of a honeycomb 
lattice with three colors in such a way that each vertex 
joins bonds of all three colors. The ground-state degener- 
acy and chiral susceptibility were obtained by Baxteoifi. 
Correlations within this set of (equally weighted) ground 
states, which determine the thermodynamics in the limit 
of low temperature, are determined by a two-component 
height model with an SU(3) symmetrjii^: one quantity 
with critical correlations is the chirality x = ±1 of colors 
around a vertes 
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FIG. 1: A portion of the kagome lattice, with the associated 
honeycomb lattice (dotted lines). For < D < J z , J xy — 0, 
and T = 0, the three spins around every small triangle include 
one with S z — +1, one with S z — — 1, and one with S z = 0. 
The bottom right shows breakup of a single bond defect into 
two vertex defects (large dots): the two defective vertices of 
the honeycomb lie on a loop containing only S z = 0, +1. 



Experimental observation of spin chirality in a kagome 
compound with approximate Heisenberg symmetry^ is 
discussed mi£. 

Now consider s = 1 spins on the sites of the kagome 
lattice, with antiferromagnetic, uniaxially anisotropic in- 
teractions J xy J z between nearest neighbors. Once 
the interactions are anisotropic, symmetry requires the 
inclusion of on-site anisotropy at the same order: 



H = 



[Jxy ■ + Sy ■ + j Z S\ ■ Si] 

(ij) 



(2) 



Note that the on-site anisotropy term D would be con- 
stant in an s = 1/2 XX Z model. Terms omitted for an 
isolated kagome layer are long-ranged (beyond nearest- 
neighbor), involve more than 2 spin operators, or break 
time-reversal symmetry, so © is appropriate when ex- 
change interactions are dominant and time-reversal is not 
explicitly broken. Dzyaloshinskii-Moriya terms absent in 
(|2"|) are allowed if the nonmagnetic lattice surrounding 
the kagome layers breaks inversion symmetry, as in the 
compound studied ini£. 

We start with the case J xy = 0, in which the model is 
effectively classical (i.e., S z on each site commutes with 
the Hamiltonian) . Rewriting the Hamiltonian as a sum 
over the small triangles in the kagome lattice gives (here 
1,2,3 label the three spins in a triangle) 
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The first term is minimized if S l z = 0, while the sec- 
ond term favors S z = for D > J z and S z = ±1 for 
D < J z . This paper concentrates on the case of weak 
positive D, but first we briefly list the other possibilities. 
For D > J z , the ground state is simply S z = every- 
where; the effective U(l) gauge theory for low-energy ex- 
7^ is added to this case has been 
For D < 0, the ground state has 
S z = ±1 everywhere and the problem reduces to the 
classical Ising kagome antiferromagnet, which is strongly 
disordered (i.e., the correlation functions in the equally 
weighted set of ground states fall off exponentially with 
distance). 

The intermediate range, < D < J Zl is classically crit- 
ical: a ground state has on every triangle one spin with 
S z = 1, one with S z = 0, and one with S z = —1, and 
this condition is an exact statement of the three-color 
problem for the bonds of the honeycomb lattice. Even 
though the Hamiltonian H does not have the additional 
Z3 color symmetry of the three-state Potts antiferromag- 
net or the classical XY antiferromagnet, the ground state 
family of states does have this symmetry. Note that the 
additional symmetry is Z3 rather than the permutation 
group of three colors S3 because there is always a time- 
reversal symmetry in the Hamiltonian: the full permu- 
tation group of the three colors is made of this twofold 
symmetry that interchanges S = +1 and S z = —1, plus 
three chirality-preserving rotations of the colors (isomor- 
phic to Z3). The icelike T — > entropy per triangle of 
the kagome lattice is half the entropy per hexagon of the 
three-coloring problem, 5a = 0.189557/caiS. 

The key difference between this s = 1 realization of the 
three-color model and the standard classical XY realiza- 
tion is the existence of an energy gap above the ground- 
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state manifold in the s = 1 model. The existence of a 
critical line (a range of values of D over which the model 
is critical) is simply a consequence of this gap: the set of 
ground states is exactly the same for all values of D, but 
the energy gap to defects varies with D as now calculated. 
Clearly the energy gap to these defects collapses at the 
two endpoints: below we obtain the correlation length at 
low temperatures induced by these defects, for which it 
is necessary to understand the nature of the long-ranged 
interaction between defects. 

At nonzero temperature above this critical line, the 
free energy and correlation length can be determined us- 
ing results on defects in the three-color model, although 
there are several signatures of the breaking of Z3 color 
symmetry once states other than the ground states are 
considered. Even though the lowest energy defect in a fi- 
nite system consists of flipping an S z = site to S z = ±1, 
so the energy gap is D, the free energy correction at finite 
temperature does not go simply as exp(— D/kT) because 
this local bond defect can break up into two vertex de- 
fects (Fig. 1). For example, if we flip S z = site to 
S z = 1 site, then by switching S z = 1 and S z = sites, 
this bond defect can fractionalize into two (1, 1,-1) ver- 
tex defects. This pair of vertex defects always lie in a loop 
made of S z = and S z — 1 sites, and are neither tightly 
bound nor completely free, but have a power-law entropic 
interaction obtained by the Coulomb gas method^ and 
verified in numerical transfer-matrix studies^. 

For the model the singular part of the free en- 
ergy per site and correlation length are related to the 
scaling dimension of lowest energy defects. The scaling 
dimension x of the defect can be calculated from the cor- 
responding contour loop model, which gives x — 
Following 1 ^^, the scaling of free energy can be calcu- 
lated in terms of scaling dimension of defects, 



f s (T) ~ \Y{m)Y(-m)Y- 



(4) 



where Y(m) and Y"(— m) are the fugacities of a pair of op- 
positely charged defects. Hyperscaling predicts that the 
singular part of the free energy per volume is determined 
by the inverse correlation volume: 
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Jxy/Jz 



f a (T) ~ £(ZT 
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FIG. 2: Phase diagram of spin-1 XX Z kagome antiferromag- 
net for fixed D £ (0, J z ). The critical point C separates 
two distinct ordered phases Pi and P2 with sixfold symme- 
try breaking. The correlation length along the dotted line 
diverges as T — > according to equation J|)J. 



with system size: in the thermodynamic limit, this pref- 
actor gives rise to the different scaling in (JSJ). 

Recent interest in the classical three-color model has 
focused on its unusually slow dynamics at low temper- 
ature, which invalidate local Monte Carlo methodsi2*2ii: 
the s = 1 anisotropic kagome antifcrromagnet may pro- 
vide an experimental system in which the dynamical 
phenomena predicted for this classical model can be 
observed. The rest of this paper concentrates on the 
quantum model: nonzero J xy generates two different 
quantum ordered phases separated by the classical line 
(Jxy — 0,0 < D < J z ) with critical correlations. Unlike 
the power-law correlations on the critical line, the or- 
der parameters in these phases are predicted to survive 
to T > 0. The following section starts with a heuristic 
derivation of the ordered phases, then develop a lattice 
U(l) gauge theory to justify some assumptions and con- 
nect to other models. 



III. PLAQUETTE-ORDERED PHASES NEAR 
THE CRITICAL LINE 



Hence the free energy and correlation length behave at 
low temperature as 



fs(T) 



-2D/3kT 



^D/3kT 



(6) 



(These formulas describe the exponential part of the be- 
havior: it is possible that there are power-law depen- 
dences in the prefactors.) The chirality correlator (JJJ 
will show power-law behavior for r = |r*i — r*2 1 *C £, and 
then decrease as exp(— r/£) once r ^> £. A small numer- 
ical transfer-matrix study (up to 9 bonds per unit cell) 
confirms qualitatively the prediction that f s for a small 
system of size L < £(T) goes as e~ D ^ kT when 
with a temperature-dependent prefactor that increases 



This section starts by finding states near the critical 
line using a standard perturbation-theoretic approach 
within the degenerate manifold of ground states found 
above. The first splitting of the ground-state degeneracy 
of the three-color model occurs at third order in pertur- 
bation theory in J xy / J z : the effective Hamiltonian is 

H eff = -t^2(SlS 2 _ SlSiS^St +h.c.) (7) 
o 

where t = (3J xy )/(2J z ), and the numbers 1 to 6 refer to 
consecutive spins around a hexagon of the kagome lattice; 
spins live on bonds of the hexagons in Fig. 1. Similar 
loop terms occur in many models when an effective H 
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is derived for states in the constrained space. A similar 
perturbation theory has been carried out for the XXZ 
antiferromagnet on the triangular lattice^. 

H e ff will induce resonances when a hexagon is "flip- 
pable" , i.e., will superpose the two configurations S\ — 
(f , 0, f , 0, 1 , 0) and (0, 1, 0, f , 0, 1 ) around a hexagon, or 
superpose (-1,0,-1,0,-1,0) and (0,-1,0,-1,0,-1). 
Note that the third type of alternating hexagon, e.g. 
(—1,1,-1,1,-1), is not flippable, so that at this lowest 
nontrivial order of perturbation theory the Z3 symme- 
try of the classical critical line is broken. The existence 
of two kinds of hexagons that can resonate will lead to 
an additional symmetry breaking in the ordered phases 
we find, compared to predictions for the s — 1/2 XXZ 
models and the honeycomb quantum dimer model&Si 

For nonzero t, we start by assuming that the system 
will favor a state with as many flippable hexagons as 
possible. This leads to the so-called \/3 x \/3 state, in 
which 2/3 of the hexagons are flippable, as a classical 
starting point. The hexagons in the honeycomb lattice 
form a triangular lattice with three sublattices: in one 
V3 x \/3 state, for example, sublattice A has spins (1,- 
1,1,-1,1,-1), sublattice B has spins (1,0,1,0,1,0), and C 
has (-1,0,-1,0,-1,0). Now all hexagons of sublattices B 
and C are flippable. To minimize the energy Q start- 
ing from this state, one expects resonance in the B and 
C sublattices, but both cannot simultaneously resonate 
without violating the three-color constraint. Consider a 
variational wavcfunction in which only sublattice B res- 
onates. The spins with S z = — 1 are fixed and do not lie 
on the hexagons of sublattice B: call these sites B. The 
other spins resonate: the lowest-energy superposition is 




O = 



± 



o = 



Oeb 



V2 




(8) 



Here the sign in the superposition is the same for all 
hexagons, and for lowest energy has the same sign as t 
and J xy : we label the state with t < as PI and that 
with t > as P2. This state has (H e ff) = — \t\/3 per 
hexagon of the whole lattice, since only 1/3 of hexagons 
resonate. It is an eigenstate of H e f f if the Hilbert space is 
restricted to states satisfying the three-color constraint. 

These states correspond to sixfold symmetry breaking: 
hexagons on one of three sublattices, containing either 
S z = +1,0 or S z — 0,-1, resonate, and the other two 
sublattices are symmetric. This combined sublattice and 
time-reversal symmetry breaking is expected to survive 
to finite temperature, with a transition temperature T c 
of order \t\ (Fig. 2). At finite temperature or by includ- 
ing additional terms in the model, it may be possible to 
separate this sixfold symmetry breaking into two transi- 
tions. The zero-temperature transition between PI and 
P2 is unusual in that it is a first-order transition by most 
definitions (the state changes discontinuously at t — 0), 



but with algebraic correlations at the transition point de- 
termined by the SU (3) theory of the three-color model: 
effectively the kagome lattice is fine-tuned so that this 
scenario occurs. For integer spin s > 2, the three-color 
problem will again describe the J xy — model if suffi- 
cently strong S z anisotropy is added to the Hamiltonian. 

A theory of two compact U(l) gauge fields on the 
bonds of the honeycomb lattice gives insight into how 
the ordered phases and critical line of this model con- 
nect to recent work on 2D and 3D frustrated magnetism. 
It is possible to use a single U(l) gauge field, follow- 
ing the compact QED description of the quantum dimer 
modeiSi: an integer-valued variable E % — S\ is assigned 
to each bond i of the honeycomb lattice to represent spin, 
with a gauge constraint that the sum of integers around 
a vertex is 0, and an energetic constraint that restricts 
E' L = —1,0,1. However, the theory now given in terms 
of two gauge fields better captures the symmetries of the 
three-color line at t = 0: the two gauge fields are related 
to the two components of the height-model description. 

First let two scalar fields E l a , a — 1,2, be assigned to 
each bond i: the values for the three spin possibilities are 



(1,0) ifS*=0 
(Ei,E*)={ (-1/2.V3/2) if 5^ = 1 

(-l/2,-V3/2) if S* = — 1. 



(9) 



These three possible states can be selected by allow- 
ing (E\,E2) to range over the sites of a triangular lat- 
tice centered on the origin with an additional energy 
\ [( E i) 2 + ( E 2?] , * 0. The Hilbert space is con- 
strained to have the sum of E fields vanish for the bonds 
around a vertex, creating two U{1) gauge symmetries. 

Next, note that a 2D unit vector can be assigned 
parallel or antiparallel to each bond i so that vertices of 
one sublattice of the honeycomb have 3 incoming bonds, 
while those of the other sublattice have 3 outgoing bonds. 
Now define two-component vector fields on bonds: E a = 
E a hi, a = 1,2. Then the three-color constraint is Gauss's 
law with no background charges: V • E Q = 0. Two height 
fields h a on the faces are defined so that scalar E a on a 
bond is equal to h a — where (r, I) indicate right and 
left w.r.t. the orientation on bonds: E a = (z x n) • Vh a , 
which guarantees Gauss's law around each vertex. This 
definition of height variables corresponds to the standard 
definition in the classical 3-color model 

The last step is to represent the ring-exchange terms, 
which break the Z3 symmetry. Define conjugate opera- 
tors on each bond (A\, A\) with commutation relations 



[Ai,E^]=tS jk S ab . 



(10) 



Then T£ = exp(iA J a ) acts as a raising operator: it in- 
creases the quantum number E J a by 1. This enables 
a compact representation of the ring-exchange terms 
proportional to t: on bond j, exp(i^/i 1 '') will raise 
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Si = to S 3 Z = 1 if = (-3/2, a/3/2). Similarly, if 

«( 2 ) = (-3/2,-^/3/2) then exp(iA^i 2) ) takes 5| = to 
= — 1. Defining vector fields = A J a hj, the ring- 
exchange terms around hexagons become 

H' = -It (cos(V x AJ { a 1] ) + cos(V x A a l 1 ^)) . (11) 

Here, as usual in gauge theories of lattice spin models, the 
meaning of cos(V x A) is that one takes the lattice circu- 
lation around a plaquette: for Vj an clockwise assignment 
of unit vectors along the bonds around a hexagon, 

cos(V x A) = cos(Ai - A 2 + A 3 - A 4 + A 5 - A 6 ) 

= cos ( J2 v, • A, j . (12) 

Note that the two ring-exchange terms only act in one 
of two cases: either the spins around the hexagon alter- 
nate between S z = Q and S z = 1, or between S z = 
and S z = — 1. The sixfold symmetry breaking in the 
plaquette-ordered phases results from the absence of the 
possible ring-exchange term with = (0, \/3). Adding 
this term restores a Z3 color symmetry, with likely nine- 
fold symmetry breaking in the ordered phases. The simi- 
larity of this gauge model to compact QED supports the 
conclusion that there is a gapped ordered state for t ^ 0. 

A quantum height model can be obtained from this 
gauge theory: using the definition of height vari- 
ables given above, the ring exchange term becomes 

Yli=i ~ 2i cos(7T a ia ), where 7r a is the operator conjugate 
to height. The quantum height Hamiltonian can be used 
to study a simple variational wave function that gives 
the same ordered state as (JSJ ■ A generalization that can 
be studied using this technique is to an s — 1/2 XX Z 
kagome bilayer in a magnetic field, with coupling between 
the layers that induces a three-color point: this system 
has columnar order rather than plaquette order. 



IV. CONCLUSIONS 

We have used the proximity to a nontrivial classically 
solvable line to study part of the phase diagram of the 
s = 1 antiferromagnet on the kagome lattice, including 
phases where J xy ^ and quantum Monte Carlo stud- 
ies are impractical. The solvable line at zero tempera- 
ture has an emergent Z3 symmetry of the three colors 
that does not exist in the full spectrum. A promising 
way to extend this approach closer to the isotropic point 
(J xy — J Z ,D = 0) is via numerical series expansion in 
Jxy/Jz from the critical line, to map out the extent of 
the plaquette-ordered phase. Finding other cases where 
frustrated magnets have solvable critical points for finite 
s similar to those in the kagome and pyrochlore antifer- 
romagnets at s = 00 is an important direction for future 
research. 

The unusual phase diagram of the spin-1 antiferromag- 
net on the kagome lattice should be verifiable by exper- 
imental neutron scattering studies of such compounds: 
in particular, confirmation of the plaquette ordered state 
should be feasible if it survives for a significant range of 
values of J xy /J z . The special low-temperature physics of 
the three-color model at J xy — also can be discerned by 
sensitive thermodynamic measurements in principle, but 
away from the special three-color point, the zero-field 
thermodynamic signature of the phase transition into 
the plaquette-ordered phase will be similar to standard 
second-order phase transitions at finite temperature. 
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